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The half-skyrmions that appear in dense baryonic matter when skyrmions are put on crystals 
modify drastically hadron properties in dense medium and affect strongly the nuclear tensor forces, 
thereby influencing the equation of state (EoS) of dense nuclear and asymmetric nuclear matter. The 
matter comprised of half skyrmions has vanishing quark condensate but non-vanishing pion decay 
constant and could be interpreted as a hadronic dual of strong-coupled quark matter. We infer 
from this observation a set of new scaling laws - called "new-BR" - for the parameters in nuclear 
effective field theory controlled by renormalization-group flow. They are subjected to the EoS of 
symmetric and asymmetric nuclear matter, and are then applied to nuclear symmetry energies and 
properties of compact stars. The changeover from the skyrmion matter to a half-skyrmion matter 
that takes place after the cross-over density ni/2 makes the EoS stiffer and leads to a compact star 
as massive as ~ 2AMq. Cross-over density in the range 1.5no ^ 2.0no has been employed, 

and the differences between the EoSs before and after this density and the dependence of our results 
on ni/2 are discussed. The novel structure of the EoS obtained with the new-BR scaling is relevant 
for neutron-rich nuclei and compact star matter and could be studied in RIB (rare isotope beam) 
machines. 

PACS numbers: 



1. INTRODUCTION 

The topological soliton called skyrmion fT] has turned 
out to be exceedingly pervasive in a variety of space- 
time dimensions ranging from 3 to 5 in many areas of 
physics and has been beautifully observed in such sys- 
tems as quantum Hall or cold atoms and more recently 
in a monoatomic magnetic film (see e.g., Q). In con- 
trast, the situation with its role in nuclear physics has 
been much less clear and with rather limited success. In 
this note, we make an attempt to uncover the power, 
hitherto unexploited, of skyrmion in strong interaction 
physics focusing on nuclear and dense matter. In contrast 
to condensed matter, the effect of skyrmion structure in 
strong interactions turns out to be indirect and hence 
much less transparent. In this work we show with simple 
but plausible arguments that the skyrmion picture can 
indeed make a clear-cut prediction on the properties of 
compact stars that has not been made thus far by other 
approaches. 

The arguments made in formulating the theoretical 
framework are neither rigorous nor completely unam- 
biguous. We rely on what Nature indicates at normal 
densities and then extrapolate to high densities using 
hidden local symmetry (HLS) structure with well-defined 
degrees of freedom . 

The starting point of our work is that when a large 
number of skyrmions as baryons are put on an FCC 
(face-centered-cubic) crystal to simulate dense matter, 
the skyrmion matter undergoes a transition to a mat- 
ter consisting of half-skyrmions Q in CC configuration 
at a density that we shall denote as This density 



is difficult to pin down precisely but it is more or less 
independent of the mass of the dilaton scalar, the only 
low-energy degree of freedom that is not well-known in 
free space. It has been estimated to lie typically at be- 
tween 1.3 and 2 times the normal nuclear matter density 
no ^5] . 

The half-skyrmion phase, made up of fractionized 
baryon numbers, is characterized by the quark conden- 
sate (qq) that vanishes on the average in the unit cell 
with, however, chiral symmetry still broken, so the pion 
is present. Note that it will have an inhomogeneous spa- 
tial distribution of baryon density. There is no obvious 
order parameter for the "transition" although there can 
be higher-dimension field operators that could be iden- 
tified at quantum level. What can distinguish the two 
"phases" are the different degrees of freedom with differ- 
ent topological charges. 

Among the predictions made so far with the half- 
skyrmion phase, one of the phenomenologically interest- 
ing predictions was the propitious drop in mass of the 
anti-kaon propagating in medium at the density ni/2 
which can have a nontrivial consequence on kaon conden- 
sation in compact-star matter Another prediction - 
which is the main object of this article and has not been 
made by other approaches - was that the presence of 
ni/2 strongly modifies the tensor forces in nuclear interac- 
tions and in particular the symmetry energy at densities 
n > ni/2 0)[a|- In this note, we confront our predictions 
with nature by translating the (semi-) classical results of 
0, 0| into the parameters of an effective Lagrangian hav- 
ing chiral symmetry and conformal symmetry, and then 
do a quantum-EFT calculation for nuclear matter and 
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compact-star matter using a renormalization-group (RG) 
based formalism 0. For the range nxji = (1-5 — 2.0)rio 
considered in 0, we have applied our formalism to neu- 
tron star calculations and a comparison of our results 
with the recently discovered two-solar-mass neutron star 
[loj will be discussed. An interesting result of the calcula- 
tion is that the skyrmion-half-skyrmion crossover makes 
the EoS stiffer at the crossover density riij^ and beyond, 
thereby leading to more massive stars. 



2. NEW-BR IN THE HALF-SKYRMION PHASE 

What plays a key role in our development is the nuclear 
symmetry energy computed in [3] in dense skyrmion mat- 
ter. There it was found that the symmetry energy Esy^n 
figuring in the energy per particle of asymmetric nuclear 
matter at density n in the form 

E(n, a) = E{n, a = 0) + Esy,n{n)a^ + 0{a^) (2.1) 

where a — {rin — rip) / (n„ + Up) with nn{np) the number 
density of neutrons (protons) is given by 



E. 



sym 



1 

8A7 



(2.2) 



where A/ is the isospin moment of inertia obtained by ro- 
tational quantization of the multi-skyrmion system which 
is given by an integral over the unit cell of a certain 
combination of the skyrmion configuration. We should 
understand that this is a classical potential energy con- 
tribution and contains no kinetic energy term. In what 
follows, we shall take into account quantum corrections 
arising from nuclear correlations. For the moment we fo- 
cus on (|2.2p . A striking feature of (|2.2p discovered in 
is a cusp at ni/2 of the symmetry energy which decreases 
from no to ni/2 and then increases for n > ni^2- Now 
given the classical nature of (|2.2p and the neglect of the 
kinetic energy term, one cannot expect this feature to 
appear unscathed in experiments. In order to confront 
nature, one has to go beyond the classical approximation 
of the skyrmion crystal. Unfortunately we are unable 
to systematically make quantum corrections within the 
skyrmion crystal approach. What we shall instead do 
is to "translate" the classical result of into the frame- 
work of an effective field theory treated at mean-field of a 
HLS Lagrangian ll| that contains all relevant degrees of 
freedom at the energy scale involved, i.e., baryons, pious 
and vector mesons. In addition, a dilaton scalar denoted 
X is introduced to account for the spontaneously broken 
conformal symmetry as precisely defined in [12:]. The 
work in ^ uses the nonlinear sigma model, involving only 
pious and baryons (emerging as skyrmions) . However the 
nonlinear sigma model Lagrangian is gauge-equivalent to 
the HLS Lagrangian, hence can capture the physics of 
vector mesons as was proposed in [l^. Our strategy 
which is consistent with the spirit of the renormalization 
group is then to do effective field theory calculation with 



this Lagrangian, with the parameters of which "running" 
with the intrinsic medium dependence as formulated in 
p^ . We will refer to this medium dependence as "BR 
scaling." 

We will now describe how the cusp structure in (|2.2p 
can be reproduced by an effective Lagrangian in mean 
field. 

Up to the density ni/2, our effective Lagrangian will 
carry the parameters scaling as introduced in [l3|. Let 
us call it "old-BR." They are of the form 



m*y/mv « m^/mjv « (x)7(x) « .K/U = (2.3) 



and 



9*/g « 1, 



(2.4) 



where the asterisk represents density dependence, 
is the pion decay constant, the subscripts N and V 
stand, respectively, for the nucleon and the vector mesons 
V = p,oj and g is the hidden gauge coupling constant . 
We will generally assume that the flavor U{2) symmetry 
applies to (p,w) in baryonic matter up to ni/2 as it does 
in matter-free space ."^^ However at n ^ '^■1/21 the frac- 
tionization of the skyrmions produces a change in the 
intrinsic scaling as [7| 
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V{n) (2.5) 



where y{n) is an order 1 constant that is more or less 
density-independent as explained below. The scaling 
is unknown except (perhaps) very near chiral transition. 
It needs not scale in the same way as <&/ does as explained 
below. Very near chiral transition, however, the HLS 
theory has, whether viable or not, a definite prediction 
thanks to the "vector manifestation fixed point (VM)" at 
which the matching of both the vector and axial- vector 
correlators gives [11] 



'/mv ^ g* /g ^ {qq)* /{qq), n 



(2.6) 



We win call (1131) with (E^]) "new-BR." 
There are two points to note here: 



• One is that at rti/2 the scaling parameter changes 
from the pion decay constant scaling as (/^//jr)^ ~ 
{qq)* I {qq) to the hidden local symmetry coupling 
constant g scaling linearly as in (12. 6p . This 
changeover was already observed in il4] from phe- 
nomenology. Since g is directly connected, via 
renormalization group flow, to the quark conden- 
sate which is the bona-fide order parameter of chi- 
ral symmetry, near the VM fixed point, it is the 



As will be elaborated below, we have found that Nature does not 
seem to agree with this symmetry for the coupling constant at 
high density. An tj-NN coupling of the form of 112.61 1 suggested 
by p-NN coupling would lead to an EoS which becomes much 
too soft in both symmetric and neutron-rich matter for density 
n > 3no. 
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vector meson mass cx g that carries information 
on chiral symmetry, not the pion decay constant. 
This changeover of the scaling from /* to g* that 
accounts for that $/ and $// need not have the 
same behavior in density is reflected in the half- 
skyrmion phase in that the pion decay constant 
drops only slowly in contrast to the quark conden- 
sate which drops to zero at ni/2- It is important to 
note that the hidden gauge coupling g scales in the 
same way as the vector meson mass does at high 
density whereas at low density up to ni/2, the gauge 
coupling stays unsealing. This difference will turn 
out to have a drastic effect on the p tensor force for 
n > ni/2 

• The second point, also connected to the slowly 
dropping pion decay constant, is that the nucleon 
mass scales little beyond 71.1/2, remaining non-zero 
at the chiral transition. This resembles - and 
we believe is related to - the nucleon mass in 
the parity-doublet nucleon model where there is 
a rather large chirally invariant mass ttiq that re- 
mains at the transition [l5| . In our application 
to be given below, we will consider mo ~ (0.7 — 
0.8)mN. 

That the new-BR affects the nuclear tensor forces 
across the density ni/2 was explained in 0, Q. So we 
will skip the details and briefly summarize only the main 
features that we will need below. 

If one takes the nucleon to be heavy while other 
hadrons, i.e., mesons, are light, then one can take the 
nonrelativistic approximation for the nucleons and write 
the effective tensor forces in medium in the usual form 
with the parameters of the Lagrangian carrying the in- 
trinsic density dependence a la BR scaling. The two ten- 
sor forces contributed by the pion exchange and the p 
exchange are given in the standard form with the masses 
and coupling constants replaced by the starred quanti- 
ties. Dividing by the vacuum quantity (Cm)^ = i^^^)"^ 
and writing = rn^jr, we have the in-medium tt and 
p tensor forces in the form 



Vlj{r)/{CMf = SmTi ■ T2S,2iR*Mfm*jY{xl,) (2.7) 



with 
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and RIj 



where M = n, p, 

Imm/Imn- The crucial (very well-known) thing to note 
is that the two forces come with an opposite sign. 

As argued in ?| and generally adopted by others, the 
pion tensor force can be taken unsealing in all relevant 
density range. In fact, one can verify explicitly that using 
suitably scaling parameters for all parameters that enter 
in the pion tensor force, such as etc., gives results that 
are close to those obtained by taking all the parameters 
unsealed [ij. Thus K « 1 and m* « in Eq. ^TT^. 



As for the p tensor, what remains to be determined is 
the scaling of i?*. It follows straightforwardly from 
and (E31) that 



g*mymN 
gmvin*^ 
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for ^ n <^ ni/2 



(2.9) 



and 



R*n 



gmvm]^ 



for 



il/2 



< n <^nc 



(2.10) 



where ric is the putative chiral transition density. Since 
g*/g ^ ^Yj in Region II, we have (|2.10p with i?* scaling 
as ^ i^Jj)^ for n ^ ni/2 and this makes a big change 
in the behavior of the net tensor force. Up to density 
n « no, the scaling $/ makes the total tensor strength 
weakened at increasing density because the increased p 
tensor eats into the pion tensor d, [l^. Recently this 
mechanism has been shown to explain the long standing 
problem of the carbon 14 dating [l7l | which in turn de- 
termines how $7 scales up to no."^"^ If the scaling $7 
continued beyond uq, then it would make the net tensor 
attraction vanish at n ~ 2no . Now with the new scal- 
ing, this behavior no longer holds. The simple prediction 
is that the net tensor-force strength will cease to drop 
at ni/2. Just to have an idea of what this does, take 
$77 $7. A simple estimation shows that when den- 
sity reaches n (2 — 3)rio, the p tensor becomes totally 
negligible. What remains is only the pion tensor. When 
this happens, tt^'s could condense in a crystalline form 
as suggested in 0, II] ■ 

We now argue that an effective field theory at mean 
field with the tensor force that follows from the new-BR 
can reproduce the cusp structure in the symmetry energy 
(|2.2|) seen in the skyrmion-crystal calculation. This can 
be seen from the fact that the symmetry energy is domi- 
nated by the tensor force E^l • ^ simple formula that 
captures the essential physics of the tensor forces is that 
of Brown and Machleidt [l^ that we rewrite including 
the new-BR, 



E, 



(2.11) 



where C is a known constant, E is the average energy 
appropriate for the tensor force, « 200 MeV, Vr is the 
radial part of the tensor force that includes the effect of 
the new-BR. There is again the kinetic term which we 



It has been suggested that this could also be explained by cer- 
tain short-range three-body forces I20II . This of course does not 
mean that three-body forces are an alternative to the scaling 
mechanism. More on this point in Discussions section. 
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will ignore as before. In the form of (|2.11l) . the cusp 
structure then follows immediately from the discussions 
given above for the behavior of the tensor forces across 
ni/2, i.e., the decrease of the net tensor force strength 
from no to ni/2 and its increase after rti/2 with the pion 
tensor taking over the strength. We take this as a support 
for the skyrmion crystal - mean-field EFT transcription. 

We will see below how nuclear correlations that go 
beyond the mean-field approximation modify this cusp 
structure. It actually smoothes it without completely 
eliminating it. The effective field theory anchored on the 
Lagrangian endowed with the new-BR is applied first to 
nuclear matter and then to compact-star matter address- 
ing the issue of the maximum neutron-star mass vs. ra- 
dius. 



3. NUCLEAR EQUATION OF STATE 

So far we have been discussing qualitative features im- 
pacted by the new-BR. We now confront quantitatively 
the scaling relations ((2?3| . ((231) . ((2J| and (p?TO| with the 
properties of symmetric as well as asymmetric nuclear 
matter. We incorporate the new-BR in the nuclear ef- 
fective field theory in which the RG-implemented Viow-k 
plays a key role. More specifically, we apply the new-BR 
scalings discussed in Section [5] to nuclear matter, both 
symmetric and asymmetric, and to the nuclear symme- 
try energy. The special features in the new-BR are the 
scaling of the nucleon mass y{n) and that of the vector 
coupling g. As we shall elaborate in the discussion sec- 
tion, these features are thought to be closely connected 
to how (most of) the nucleon mass is generated in the 
strong interactions. 

Before continuing, let us concisely recapitulate how 
the new-BR enters into an RG-implemented EFT. As 
argued in it involves two decimations in the RG 

sense. Starting with an effective chiral Lagrangian, one 
first decimates in matter-free space from the chiral scale 
A^ ^ Airf-n ^ 1 GeV down to the first decimation scale 
A ^ 3 fm^^. What results is the Viow-k that is used in 
our calculation. Then doing many-body calculations for 
nuclear systems with the parameters of Viow-k running 
a la new-BR amounts to doing the second decimation. 
In fact this second decimation is equivalent to doing a 
Landau Fermi-liquid theory calculation as formulated in 
14| . In doing this, we are ignoring 3-body and higher- 
body forces. One should however recognize that part of 
many-body force effects are embedded in the new-BR. 
One can think of this as a sort of duality between the 
two as will be elaborated later. 

To suitably take into account the features mentioned 
above into a high-order effective field theory calculation, 
we shall carry out our calculations using the realistic 
BonnS potential [2l[; this potential is an extension of 

111 with the 



the one-boson-exchange BonnA potential 
provision that the nucleon and meson mass as well as 
the vector coupling be scaled a la new-BR. As discussed 



in Section [2l we employ the following two-region scal- 
ings characterized by the transition densities ^tnd ric 
(respectively for the skyrmion-half-skyrmion and chiral 
transitions). For density < n < ni/2 (Region-I), we 



and 



= ^ = 1>/(n);<I>/(n) 
rriM TTiN 



R* = 1. 



l + ci^ 

^ no 



(3.1) 



(3.2) 



In the above M ~ {V, S) and N stand respectively for 
meson (both vector and scalar) and nucleon. For density 
ni/2 < n <nc (Region-II), we use 



niM 

for mesons and 



<I>r,(n);$-(n) = 



1 + 

no 



^ = $-(n)=y(n) 
ruN 

for nucleons. We use the R* scaling in II as 



9 



(3.3) 



(3.4) 



(3.5) 



The above scaling functions and $// are in general not 
continuous at the boundary density 71-1/2. This disconti- 
nuity may be a mere artifact of the simplification we are 
adopting. In the present work, as to be discussed later, 
we shall choose the parameters contained in them so that 
these two functions are nearly continuous (to avoid dras- 
tic discontinuity) at ni/2- In addition, we shall employ 
two Fermi-Dirac functions to smoothly join the scaling 
functions and so that the resulting scaling func- 
tion $ is ensured to be continuous at the boundary. A 
similar procedure will also be employed for the R* scal- 
ings in the two regions. To illustrate, the smoothed scal- 
ing function $ is constructed as 



with 



$ = ^^<(ni/2)$/-f ^^>(ni/2)$77, (3.6) 

F<{ny2) = [l + e("-"V2)/5]-i^ 

F>{ny2) = [l + e("V2-«)/5]-i (3.7) 



We must stress that except for low density < no (and possi- 
bly high density near the chiral transition point in the chiral 
limit as predicted in HLS), the precise form of the scaling is 
not known, so what we take should be understood as more of 
a convenient parametrization guided, whenever feasible, by phe- 
nomenology. Furthermore there is nothing that suggests that the 
scaling should be identical for all mesons. 



where ^ is a smoothness parameter. In the present work, 
we shall use S/uq ~ 0.05-0.10. It turns out that withm 
this range our results are satisfactorily stable with respect 
to 5. 

We have adopted the following procedure for choosing 
the parameters of the above scaling functions. First we 
require the parameters in Region I so that they satisfac- 
torily reproduce the empirical nuclear matter saturation 
properties (saturation density ~ 0.16 fm^'^ and av- 
erage energy per nucleon £^o/^ ~ — 16 MeV at satura- 
tion). The choice for the parameters in Region II will 
be addressed later. We shall calculate no and E'nM us- 
ing a low- momentum ring-diagram approach [23l - l26l | , 
where the pphh ring diagrams are summed to all orders 
within a model space of decimation scale A. Few low- 
order (1st-, fourth- and eighth-order) such diagrams are 
displayed in Fig. 1. Note that each vertex of the diagrams 
is a low-momentum interaction Viow-k which is obtained 
from a realistic NN potential Vnn using a renormaliza- 
tion group approach where the momentum components 
beyond a decimation scale A are integrated out |27H30| . 

More precisely Viow-k is given by the following T- 
matrix equivalence equations: 



T{k',k,k^) = VNN{k',k) 



TT 



FNN(fc',<z)r(g,fc,fc2) 



Tiovj—k{k , fc, k 

+ -r 

TT 



') ='l^ow-k(fc',fc) 

V^ow-k(fc^ g)riow-k(g, fc, fc^) 
fc2 



q'^dq, (3.. 



q^dq, (3.9) 



T(fc', fc, fc2) = Tion,~k{k', k, fc2); (fc', k) < A. (3.10) 

In the present work the above Vnn is chosen to be the 
realistic BonnS [21*] NN interaction. (The new-BR scal- 
ings we have established above enter into the meson pa- 
rameters as well as the nucleon mass of this potential 
with the varying density.) V denotes principal- value in- 
tegration and the intermediate state momentum q is in- 
tegrated from to cxD for the whole-space T and from 
to A for Tiow-k- Because we shall calculate the nuclear 
symmetry energy Esym{n) up to n ~ 5no, we shall use 
A = 3 fm~-^ p. The above Viow-k preserves the low- 
energy phase shifts (up to energy A^) and the deuteron 
binding energy of Vnn- (For example, the deuteron bind- 
ing energy given by Viow-k of A = 2.0 and 3.0 fm~^ are 
both -2.226 MeV.) Since Viow-k is obtained by integrat- 
ing out the high- momentum components of Vnn, it is 
a smooth 'tamed' potential which is suitable for being 
used directly in many-body calculations. The familiar 
HF approximation for nuclear matter corresponds to the 
inclusion of only the first-order diagram (a) of the figure. 
In contrast, the pphh ring diagrams such as those shown 
in Fig. [T] are included to all orders in our nuclear matter 
calculations. 



0--0 




(a) 



(b) 



FIG. 1: Diagrams included in the all-order pphh ring-diagram 
summation for the ground state energy of nuclear matter. 
Each dashed line represents a Viom-k vertex. 



With such ring diagrams summed to all orders [23|, |2J] , 
the ground-state energy of asymmetric nuclear matter is 
expressed as E{n,a) — E-^^'^'^{n,a) + AE{n,a) where 
jTjfree (j^gj^gtes the energy for the non-interacting sys- 
tem and AE, the energy shift due to the NN interac- 
tion, is given by the all-order sum of the pphh ring di- 
agrams as illustrated in Fig. [1] We include in general 
three types of ring diagrams, the proton-proton, neutron- 
neutron and proton-neutron ones. The proton and neu- 
tron Fermi momenta arc, respectively, kpp — {Sir'^npy/^ 
and kpn ~ (STr^nn)^/'^, where Up and n„ denote respec- 
tively the proton- and neutron-density. The asymmetric 
parameter is a = (n„ — np)/(nn + rip). With such ring 
diagrams summed to all orders, we have 



m ijkKA 

xr„;(A:Z,A)(ij|Viow-k|fcO, (3.11) 

where the transition amplitudes Y are obtaind from a 
pphh RPA equation [23l[23|. Note that A is a strength 
parameter, integrated from to 1. The above ring- 
diagram method reduces to the usual HF method if 
only the first-order ring diagram is included. In this 
case, the above energy shift becomes AE{n,a)HF = 
5 Z]"»"j(«i|Mow-k|u) where nfe=(l,0) if fc(<, >)fcFp for 
proton and nfe=:(l,0) if k{<,>)kpn for neutron. 

The above Viow-k ring-diagram framework has been 
app lied to symmetric and asymmetric nuclear matter 
[23L HJ] and to the nuclear symmetry energy . This 
framework has also been tested by applying it to dilute 
cold neutron matter in the limit that the ^Sq scatter- 
ing length of the underlying interaction approaches in- 
finity pa . [2^. This limit - which is a conformal fixed 
point - is usually referred to as the unitary limit, and 
the corresponding potentials the unitarity potentials. For 
many-body systems at this limit, the ratio ^ = Eq/Eq^'^'^ 
is expected to be a universal constant of value ^ 0.44. 
{Eq and Eq^'^'^ are, respectively, the interacting and non- 
interacting ground-state energies of the many-body sys- 
tem.) The above ring-diagram method has been used 
to calculate neutron matter using several very differ- 
ent unitarity potentials (a unitarity CDBonn potential 
obtained by tuning its meson parameters, and several 
square- well unitarity potentials) [1^ [2^. The ^ ratios 
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given by our calculations for all these different unitar- 
ity potentials are all close to 0.44, in goo d ag reement 
with the Quantum-Monte-Carlo results (see [2a] and ref- 
erences quoted therein). In fact our ring-diagram re- 
sults for ^ are significantly better than those given by 
HF and BHF (Brueckner HP) [H, HI]. It is desirable 
that the above unitary calculations have provided satis- 
factory results, supporting the reliability of our Viow-k 
ring-diagram framework for calculating the nuclear mat- 
ter EoSs. 



4. RESULTS 

We recall that the new-BR has an assumption which 
is not an immediate consequence of chiral symmetry. 
Specifically the premise of the vector manifestation as- 
sociated with hidden local symmetry states that as one 
approaches the VM fixed point, my/mv g*/g 0, 
which is not dictated by chiral symmetry alone. Here 
we will take the point of view that the vector manifesta- 
tion property is operative after the half-skyrmion onset 
density ni/2- 
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FIG. 2: New-BR EoS for symmetric nuclear matter calcu- 
lated with ni/2 = 1.5no. See text for more explanations. 

Let us first consider the EoS of symmetric nuclear mat- 
ter in the low density region (n ^?ii/2), a main purpose 
here being the choice of the c/ parameters so that the em- 
pirical saturation properties of symmetric nuclear matter 
are satisfactorily reproduced. In Fig. [2l we present our 
results for symmetric nuclear matter calculated with pa- 
rameters C7=0.130 for nucleon and p-meson, =0.121 for 
cr-meson and =0.139 for w-meson.'^'' The EoS of this 
figure gives ground-state energy per nucleon Eq/A =-15 
MeV, saturation density rt^at — 0.93no and compression 



Here we are doing some fine-tuning for a better fit but the small 
diff'erences in cj's are of course of no significant meaning. 



modulus K= 208 MeV, all in satisfactory agreement with 
the empirical values. (Here and in Fig. 2 Eq /A is the same 
as {E{n,a = 0) — rriN) of Eq.(2.1).) The above calcula- 
tion has employed ni/2 — 1.5no. As to be presented later, 
we have also carried out calculations with ni/2 = 2rio 
and the saturation properties given by them are nearly 
the same as the ni/2 = 1.5no case. A decimation scale of 
A = 3 fm^^ has been employed in the above calculation, 
and it will be used from now on. 

As discussed in [2p, l23| , the use of realistic Vn n with 
the old-BR [l^ I3TI.T32] leads to satisfactory nuclear mat- 
ter saturation properties. As indicated by Fig. 2, we see 
that the new-BR does also lead to satisfactory nuclear 
matter saturation properties, even though these two scal- 
ings are largely different. The main differences between 
them are in the scaling of the nucleon mass and the HLS 
coupling g. The nucleon mass and coupling constant g 
are both not scaled in old-BR while scaled in new-BR. 
Furthermore, the region-I and -II scalings of old-BR are 
the same, while there are important differences between 
the new-BR scalings in these two regions which we shall 
later address. 

Before proceeding to the EoS for n let us first 

discuss the scaling parameters we have employed. The 
scaling functions we have used in Region I (Eqs.(3.1-3.2)) 
are similar to those employed in the old-BR flj, [s^] and 
Ericson (ER) [s^ scalings. The ER scaling is based on 
the quark condensate relation [33 [ 



Ml 



1 



1 



(4.1) 



where S^Ar=45±7 MeV i3J]. 
hadrons in medium reads 



Then the ER scaling for 



with 



D = 



m 
m 



1 



1/3 



l + D- 



(4.2) 



Using the empirical values for 
we have £»= 0.35±0.06. Note that 



for the low-density region this relation agrees well with 
the parametrization for the old-BR [l^, Is^] 



= 1 -C— 

m uq 



(4.3) 



where C is a constant of value ^ 0.15. 

It may be noted that our new-BR in Region I (Eqs. 
(3.1-3.2)) is consistent with the above Ericson scaling in 
the n < ni/2 region if the c/ scaling parameters are cho- 
sen to have values near D/3 ~ 0.12±0.02. It is encourag- 
ing that the c/ parameters we have employed so as to give 
satisfactory nuclear matter saturation properties (Fig. [2]) 
are indeed quite close to the value of D/3 given by QCD 
theories. 

We now consider the EoS for n > nx/2- In Fig- 3 we 
present results for two choices for the half-skyrmion on- 
set densities, namely ni/2= 1-5 and 2.0no. In addition 
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FIG. 3: Comparison of the EoS for symmetric nuclear matter 
calculated with n^i^ ~ 1.5 and 2.0no. See text for more 
explanations. 



we also present the EoS (labelled (C) in the figure) ob- 
tained with the unsealed BonnS [2l| potential. As seen 
this EoS does not have satisfactory nuclear matter satu- 
ration properties; it would give saturation density much 
higher than the empirical value of ^ O.lQfm"^ as well as 
a saturation energy much lower than the empirical value 
of ~ —16 MeV. In contrast, the new-BR EoSs (A) and 
(B), respectively for ni/2= 2.0 and l.SriQ, both have sat- 
isfactory saturation properties. In calculating (A) and 
(B), we have used the same c/ parameters as listed ear- 
lier. Thus (A) and (B) are equivalent for n < 1.5no, both 
having the same saturation properties {Eo/A^-15 MeV, 
Usat = 0.93no and 208 MeV). 

Turning to the EoS in Region II (n > ni/2), we note 
from (j2.9|) and (I2.10p that the scalings for R* controlling 
the p tensor force in Region II are significantly different 
from those in Region I (see Eqs.(3.1-3.5)): Other com- 
ponents of vector-meson-exchange nuclear forces are gov- 
erned, apart from the mass scaling, by the scaling of the 
hidden gauge coupling constant g 



9* la ~ I, n< ni/2\ 



(4.4) 
(4.5) 



Were the flavor U{2) symmetry operative in Regin II, 
the scaling (|4.5p would hold for both p and w. It turns 
out, however, that if the w-nucleon coupling dropped in 
the same way as the p-nucleon coupling, nuclear systems 
would collapse in that region. We have found that the 
repulsion provided by the oj-exchange potential is sen- 
sitively dependent on the w-nucleon coupling constant, 
and a moderate dropping of this constant can drastically 
suppress the repulsion, making the system unstable at 
high densities. This signals that the coupling constant 
g* must be asymmetric in high density or as mentioned 
in Discussion, higher members of the ui mesons in the 
infinite tower in holographic QCD models that arise in 
string theory [HI could intervene in providing the neces- 



sary repulsion. In our calculation, we will take g* / g ~ 1 
for the w-exchange potential for both Regions I and II. 
We will come back to this matter at the end of the paper. 

Another difference is that the scaling of the nucleon 
mass (m^/mjv) is density dependent in Region I as is 
seen in experiments while it is equal to a constant or 
slowly varying in Region II. These differences can make 
the EoSs in these two regions significantly discontinuous 
at n — ni/2- As mentioned above, this discontinuity 
could be an artifact of our schematic treatment of the 
skyrmion-half-skyrmion transition. We have found that 
this discontinuity can be made small by suitably choosing 
the scaling parameters in Region II. We have done so, and 
the EoSs (A) and (B) presented in Fig. 3 are obtained 
with c// = cj for both (A) and (B), with y(n)=0.77 and 
— 2.0no for (A), and y(n)=0.78 and — l-5rto 
for (B). We note that both values are close to 0.80. These 
parameter choices will be referred to respectively as A- 
parameters and B-parameters. They will be used and 
tested in other calculations such as nuclear symmetry 
energies and neutron stars later on. 

It is of interest that the EoSs (A) and (B) of Fig. 3 
both exhibit a narrow 'plateau' or 'zigzag' segment near 
n ^ 2.2nQ. The occurrence of this plateau/zigzag struc- 
ture may indicate a skyrmion-half-skyrmion transforma- 
tion which is of interest for further study. '^^ This occur- 
rence could be largely due to the use of a constant (non- 
scaling) nucleon mass beyond ni/2 - that we assume - 
as is indicated in the half-skyrmion matter and predicted 
in the parity-doublet model with a large chiral-invariant 
mass Too- It may be pointed out that y(n) plays an im- 
portant role in determining the EoS of nuclear matter in 
Region II as illustrated in Fig. 3. For density n > 2.0no, 
the parameters used for the EoSs (A) and (B) there dif- 
fer only in y(n)=0.77 for (A) and = 0.78 for (B), the cn 
parameters used for them being identical. It is seen that 
this small difference has made the (A) EoS significantly 
more repulsive than (B), especially in the high density 
region. Our calculations have found that the use of a 
smaller y(n) would generally give an upward lift to the 
n > ni/2 EoS, resulting in a stiffer EoS. 

By way of heavy-ion collision experiments, there has 
been much progress in determining the nuclear symme- 
try energy Egym up to densities as high as ^ 5no KUi viN . 
Thus an application of our new-BR scaling to the calcu- 
lation of Esym would provide an important test for this 
scaling in the region with n > ni/2- (We remind the 



#^ Such a changeover is observed on crystal. It is topological and 
hence likely robust. However there is no well-defined order pa- 
rameter to signal a phase transition. The quark condensate (qq) 
vanishes on average in unit cell in the crystal description but this 
is not a bona-fide order parameter since the pion is present in the 
system, indicating chiral symmetry is not restored in the half- 
skyrmion state. It appears to be possible to construct a higher- 
field operator as an order parameter representing an emergent 
symmetry in medium. This matter will be discussed elsewhere. 
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readers that our new-BR scaling has worked rather weh 
for n < ni/2-) The nuclear symmetry energy is related to 
the asymmetric nuclear matter EoS as Eq. ()2.1|) . We have 
calculated E{n, a) for a range of a values, and from them 
we Also we use the same ring-diagram for- 

malism where the pphh ring diagrams are summed to all 
orders. Our results, labelled (A) and (B) and calculated 
respectively with the same new-BR A- and B-parameters 
mentioned earlier, are displayed in Fig. 4. It is of inter- 
est that the symmetry energies given by (A) and (B) are 
nearly identical, despite the considerable differences be- 
tween the two corresponding EoSs for symmetric nuclear 
matter shown in Fig. 3. 
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FIG. 4: Comparison of our calculated nuclear symmetry en- 
ergy (new-BR) with the empirical upper (expt-Lil) and lower 
(expt-Li2) constraints of Li et al. 139] and the empirical re- 
sults of Tsang et al. (expt-Tsang) [381 ]. See text for more 
explanations. 

Based on heavy-ion scattering experiments, Li et al. 
[36j have suggested an empirical relation 



E. 



sym{n) w 31.6(n/no)''; 7 = 0.69 



I.l, 



(4.6) 



for constraining the density dependence of the symme- 
try energy. The upper (7 = I.l) and lower (7 = 0.69) 
constraints are also plotted in the figure, labelled respec- 
tively as 'expt-Lil' and 'expt-Li2'. Also based on such 
experiments, Tsang et al. [38[ recently proposed a new 
empirical formula for the symmetry energy, namely 



Esym (^) 



Cs.k 



X 2/3 ^ 



no 



no 



(4.7) 



0.7. 



where Cs,k = 25MeV, Cs,p = 35.2MeV and 7, 
This formula is also plotted in Fig. |31 labelled as 'expt- 
Tsang'. Note that Tsang's results are very close to the 
lower constraint of Li et al. 

Returning to Fig. 4, we see that our new-BR results 
agree reasonably well with the empirical constraints on 
Esym'i for Region I our results are slightly below the em- 
pirical lower bounds while in Region II they tend to be 



closer to the upper bound. Recently Lattimer and Lam 
[39} have investigated the constraints on Egy^ and L (de- 
fined as {3/8){dEsym/du),u = n/no) at density n = no- 
The results deduced from nuclear masses, nuclear giant 
dipole resonances, astrophysics, neutron skins of the Sn 
isotopes, and other investigations exhibit wide variations, 
with Esym/MeV ranging from ~ 24 to ^ 36 and L/MeV 
from ~ —20 to ^ 100. The overlap constraints allowed 
bi 

E 



all these results are 28 ^ E,ym ^ 32 and 40 ^ L <^ 60 



Our results as given by (A) and (B) of Fig. 4 are 
E,ym/MeV ~ 26.5 and L/MeV ^ 8. They are both be- 
low the respective lower bounds of the above constraints, 
although consistent with the constraints given by the nu- 
clear masses and nuclear giant dipole resonances. 
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FIG. 5: Comparison of our calculated neutron matter EoS 
with the empirical stiff (upper box) and soft (lower box) con- 
straints of Danielewicz et al. [iol ] . 

From heavy-ion collisions, Danielewicz et al. [4^ have 
obtained constraints for the density-pressure EoS pin) of 
neutron matter up to densities ^ 4.5no. To further study 
our new-BR scaling in the high density region, we have 
calculated the neutron p{n) EoSs up to the above densi- 
ties. A comparison of our results with their constraints is 
presented in Fig. 5 where the upper and lower boxes are 
resp ectively the constraints for the stiff and soft EoSs of 
[40j . Our EoSs calculated with parameters A and B are 
denoted by 'points' and 'line' respectively. 

Before discussing them, we need to briefly describe 
some calculational details. Including the nucleon rest- 
mass energy, we first calculate the nuclear-matter energy 
density 



en 



n( 



Eojn) 
A 



uin) 



(4.8) 



with the average ground-state energy Eo (n) /A obtained 
from the ring-diagram method described earlier (sec- 
tion 3). The pressure-density EoS is then given by 



p{n) 



de{n) 
dn 



e(n). 



(4.9) 
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As we have seen in Fig. 3, our EQ{n)/A has a short seg- 
ment of 'zigzag' (plateau-hke) structure at n ~ 2.2no, 
and consequently this structure is also present in e(n). 
This feature makes it very difficult to calculate the 
derivatives of e(n), hindering the calculation of p{n). 
To circumvent this difficulty, we have employed a fit- 
ting procedure so as to have a smooth (differentialble) 
e(n). Such a smooth crossover is expected also in a hy- 
brid hadron-quark matter model mentioned below [41|. 
Li and Schulze 42| recently proposed a highly desirable 
parametrization for the nuclear-matter EoS: they have 
found that a wide range of nuclear EoSs can be fitted very 
accurately by the polytrope EoS Eo{n)/A = ain + a2n'^^ 
where ai, 02 and 03 are parameters. We have found 
that this type of polytrope EOS can fit our EOS very 
well, significantly better than the other fitting methods 
we have considered. The fitted EoS and the original one 
agree very well outside the narrow zigzag region while 
within this region the fitted one is basically a smooth 
line crossing through the plateau/zigzag segment. We 
have calculated the neutron EoS p(n) based on this fit- 
ted polytrope EoS, and the results calculated with the A- 
and B-parameters specified earlier are shown in Fig. 5. 
A similar comparison for the symmetric nuclear-matter 
EoS is displayed in Fig. 6. 
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FIG. 6: Same as Fig. 5 except that it is for symmetric nuclear 
matter. 

As shown in Figs. 5 and 6, our calculated EoSs for neu- 
tron matter and symmetric nuclear matter are both in 
satisfactory agreement with the empirical constraints of 
Danielewicz et al. [10] . Our neutron-matter EoS is some- 
what stiffer in the high-density regime than what would 
give the 2-solar mass star. There seems to be nothing 
obviously wrong with this. The model having the hybrid 
hadron-quark continuity mentioned above [41| seems to 
favor such massive stars. But it should be noted that 
at high density, there are several compelling arguments 
such as the vector manifestation and the dilaton limit of 
hidden local symmetry structure for an ultimate soften- 
ing of the repulsion as density approaches ric 0] ■ What 
is not known is at what density they become operative. 
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FIG. 7: Mass-radius trajectories of neutron stars calculated 
with new-BR scalings using ~ 2.0 (A) and 1.5no (B). 

The maximum neutron-star mass and its radius for these two 
cases are respectively (2.39 Mq, 10.9 km) and (2.38 Mq, 10.9 
km). 



In a recent neutron-star calculation using realistic NN 
potentials [23|, the effects from the 'old-BR' scaling ((23)) 
and (12. 4p applied in both I and II were found to be 
highly important for neutron stars, the maximum mass 
and its radius calculated (with, without) the inclusion 
of such effects being respectively (~ 1.8, ~ 1.2Mq) and 
(~ 8.9, ^ 7.2)km. Now the question is: What does 
the new-BR (|2.5p do to neutron stars? To address this 
question, we have calculated the properties of pure neu- 
tron stars (i.e. made of neutrons only) from the above 
neutron-matter EoSs, using the calculation procedures 
described in [IJj. In Fig. 7 we present our calculated 
neutron-star mass-radius trajectories (A) and (B), ob- 
tained respectively with the A- and B-parameters men- 
tioned earlier. The main results of these two calculations 
are nearly the same, in consistence with the close similar- 
ity between the two EoSs shown in Fig. 5. The maximum 
mass of neutron stars given by (A) and (B) are respec- 
tively 2.39 and 2.38 Mq and the corresponding radii both 
^11 km. The effects from the present new-BR scaling 
appear to be even stronger than those from the old-BR 
scaling. For example, the maximum mass obtained (with, 
without) the new-BR scaling are (^ 2A,^ 1.2Mq), the 
increase between them being significantly larger than the 
above old-BR case. 



5. DISCUSSIONS 

In this paper, we subjected the nuclear effective field 
theory anchored on RG flow, with the parameters of the 
Lagrangian sliding with density, to normal nuclear mat- 
ter and dense compact-star matter. The scaling behavior 
used here differs from the old BR scaling , in that at 
a density ni/2 > no, a topological change takes place 
from skyrmion matter to half-skyrmion matter, giving 
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rise to a modified scaling "new-BR." The changeover 
from skyrmion matter to half-skyrmion matter is char- 
acterized by a vanishing quark condensate (qq) = but 
a nonvanishing pion decay constat fjr ^ 0- Thus it is not 
a standard phase transition a la the Ginzburg-Landau- 
Wilson paradigm although two different phases are in- 
volved; it appears to involve an emergent symmetry not 
present in the fundamental theory, QCD. 

At the semi-classical approximation made in the cal- 
culation, the half-skrymions are not deconfined in con- 
trast to what happens in certain condensed matter sys- 
tems They are bound or confined, so they are not 
propagating degrees of freedom. What characterizes the 
system is that the mass of the baryon made up of two 
'bound' half-skyrmions remains more or less unsealed, 
not going to zero up to the density ric at which the 
quarks get deconfined, whereas the vector meson mass 
is expected to drop faster in half-skyrmion phase than 
in skyrmion phase. This means that the origin of the 
most, if not all, of the nucleon mass is not in the dynam- 
ical symmetry breaking of chiral symmetry, in contrast to 
the meson mass, with a substantial mass of the nucleon 
coming from a hitherto unknown source. This is simi- 
lar to what is described in the parity-doublet model of 
the nucleon [l^li^.*^ It is intriguing that the two con- 
sequences of the changeover at ni/2, namely, the drastic 
modification of the nuclear tensor force and the stiffening 
of the EoS of dense matter at seem to be hinting 

at the mechanism for the generation of ~ 99 % of the 
nucleon mass in the strong interactions. 

The salient features obtained in the RG-implemented 
effective theory approach adopted in this paper can be 
summarized as follows: 

1. Without a suitable scaling in the Lagrangian that 
figures in Viow-k (or incorporating many-body 
forces), symmetric nuclear matter cannot be stabi- 
lized at the right density and with correct binding 
energy. 

2. Our calculations have essentially two scaling pa- 
rameters: one is c/ « 0.13 for all mesons (vector 
mesons and scalar meson) and the nucleon in region 
I, and in region II we have c// = c/ for mesons and 
the vector coupling and an additional parameter 
y{n) « 0.8 for the nucleon. With these two param- 
eters, one can explain satisfactorily the saturation 
density, the binding energy and the compression 



#^ This picture is clearly at odds with the constituent quark model — 
which has a strong theoretical support from QCD in the large Nc 
limit [45l — where the ratio of the meson mass over the baryon 
mass is 2/3. Whether or not the constituent quark model is 
applicable in nuclear medium is not known, but if there were an 
mo for the quark which is not small, then it should be possible 
that the constituent quark model hold in dense medium and the 
ratio remain more or less the same. In this case, the scaling could 
be considerably different from the new-BR. 



modulus of symmetric nuclear matter as well as 
the nuclear symmetry energy, and predict the EoSs 
for symmetric and asymmetric nuclear matter at 
high density and compact-star matter. Our results 
give a good fit to all quantities that are available 
experimentally at densities up to n ^ 4no. 

3. The topology change from skyrmion to half- 
skyrmion at ni/2 changes the slope of the EoS, 
making it stiffer in the half-skyrmion phase and 
raises the maximum mass of compact stars to ^ 
2AMq. Verifying the presence and the role of the 
topology change at ni/2 should be feasible at RIB 
machines. 

To conclude, we make a few additional remarks and 
mention several caveats in our work. 

In transcribing the new-BR into nuclear EFT, we have 
not included n-body potentials in Viow-k with n > 2. 
It is known that one can obtain reasonable symmet- 
ric nuclear matter without BR scaling but with three- 
body forces [13. This suggests that the BR in the two- 
body potential and the three-body forces play roughly 
the same role in the EoS of symmetric nuclear matter 
near the saturation density. In some channels, the map- 
ping between the two could be nearly complete. They are 
not alternatives. An apt example of this dichotomy (or 
duality) is the explanation of the cl4 dating beta-decay 
by either the BR [17] (up to uq, both 'old' and 'new' are 
the same) or short-ranged three-body forces [20| . 

In principle, there should be no problem in including 
both BR and many-body forces in a way consistent with 
the tenet of chiral expansion. What one has to do in 
the presence of such n-body potentials is then to suitably 
modify the scaling properties of the Lagrangian, since di- 
rect and indirect chiral symmetry effects are compounded 
in physical quantities in a variety of different chiral ex- 
pansion schemes as illustrated in [46| . A fully consistent 
way of doing the calculation would be to have both the 
scaling and many-body potentials treated together with 
certain constraints, such as thermodynamic consistency, 
taken into account. 

We have not taken into account strangeness degrees of 
freedom - such as kaons, hyperons, strange quarks etc. - 
into the EoS for neutron-rich matter. In our formulation 
anchored on dense skyrmion matter, as described in 0, 
hyperons can enter only after kaons condense. Therefore 
the issue here is how kaon condensation can take place af- 
ter changing from skyrmion matter to half-skyrmion mat- 
ter. This problem has not been fully worked out yet. The 
only thing we can say is that the increased repulsion in 
the EoS for dense matter in the half-skyrmion phase can 
push the kaon condensation threshold to higher densi- 
ties than naively estimated in chiral perturbation theory. 
The stiffening of the EoS from soft one to hard one we 
noted at ni/2 is analogous to the stiffening at the smooth 
crossover at a density ~ (2 — 4)no from hadron to non- 
strange quark phase in the hybrid model that also yields 
the maximum star mass ~ 2.3Mq [41| . One intriguing 
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possibility is that our half-skyrmion phase is related to 
the smooth crossover regirae in (411] . In some sense we can 
consider the half-skyrmion matter as a hadronic dual to 
a strong-coupled quark matter. 

There is one caveat in our treatment which requires 
further studies. We have found that if we were to drop 
the w-nucleon couphng according to gl^jg^ ~ 9p/9p ~ 
g* / g = as one would expect if flavor U (2) symmetry 
held in Region II, the EoS would become much too soft 
above hq to be compatible with the existence of the 2- 
solar mass object observed in nature. We kept g^/g^ « 1 
while letting the uj mass scale. This does not seem to 
be consistent with hidden local symmetry which makes 
a clear-cut prediction that the repulsive core due to cj 
exchange should get suppressed at high density. 

One can think of two possibilities for addressing this 
caveat. One is that the flavor U{2) symmetry could be 
badly broken for the coupling constant in dense medium. 
The other possibility is that the short-distance repulsion 
between baryons which seems to be universal in that it is 



flavor-blind is mediated by the infinite tower of to mesons 
implied in dual gravity models in hQCD [s^. While the 
coupling to nucleons of the lowest cj meson gets sup- 
pressed as required by the vector manifestation and the 
dilaton-limit fixed point of hidden local symmetry, the 
higher-lying members of the tower could intervene to pre- 
serve the repulsion that seems to be needed. These are 
open problems to be resolved in the future work. 
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